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Abstract. According to the Gottsche conjecture (now a theorem), the degree 
]\j-d,5 q£ -j-j^g Severi variety of plane curves of degree d with 6 nodes is given by a 
polynomial in d, provided d is large enough. These "node polynomials" Ns{d) were 
determined by Vainsencher and Kleiman-Piene for S < 6 and 6 < 8, respectively. 
Building on ideas of Fomin and Mikhalkin, we develop an explicit algorithm for 
computing all node polynomials, and use it to compute Ns{d) for S < 14. Further- 
more, we improve the threshold of polynomiality and verify Gottschc's conjecture 
on the optimal threshold up to 5 < 14. We also determine the first 9 coefficients 
of Ns{d), for general S, settling and extending a 1994 conjecture of Di Francesco 
and Itzykson. 



1. Introduction and Main Results 

Node Polynomials. Counting algebraic plane curves is a very old problem. In 
1848, J. Steiner determined that the number of curves of degree d with 1 node 
through '^^'^^^^ — 1 generic points in the complex projective plane is 3((i — 1)^. 
Much effort has since been put forth towards answering the following question: 

Hovj many (possibly reducible) degree d nodal curves with 
5 nodes pass through ^^^th^) _ § generic points in ? 

The answer to this question is the Severi degree N'^'^, the degree of the corresponding 
Severi variety. In 1994, P. Di Francesco and C. Itzykson [6] conjectured that N'^'^ is 
given by a polynomial in d (assuming S is fixed and d is sufficiently large). It is not 
hard to see that, if such a polynomial exists, it has to be of degree 26. 

Recently, S. Fomin and G. Mikhalkin [71 Theorem 5.1] established the polynomi- 
ality of N"^'^ using tropical geometry and floor decompositions. More precisely, they 
showed that there exists, for every 5 > 1, a node polynomial Ns{d) which satisfies 
iV^.^ = Ns{d) for all d > 26. (The 6 = case is trivial as A^^'° = 1 for all d > 1.) 

For 6 = 1,2, 3, the polynomiality of the Severi degrees and the formulas for Ns{d) 
were determined in the 19th century. For 6 = 4,5,6, this was only achieved by 
I. Vainsencher p!3] in 1995. In 2001, S. Kleiman and R. Piene [5] settled the cases 6 = 
7, 8. Earlier, L. Gottsche [Bj conjectured a more detailed (still not entirely explicit) 
description of these polynomials for counting nodal curves on smooth projective 
algebraic surfaces. 
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Main Results. In this paper we develop, building on ideas of S. Fomin and G. 
Mikhalkin [7], an explicit algorithm (see Algorithm [T]) for computing the node poly- 
nomials Ns{d) for arbitrary 6. This algorithm is used to calculate Ns{d) for all 
6 < 14. 

Theorem 1.1. The node polynomials Ns{d), for S < 14, are as listed in Appendix\A\ 

P. Di Francesco and C. Itzykson [6] conjectured the first seven terms of the node 
polynomial Ns{d), for arbitrary 5. We confirm and extend their assertion. The first 
two terms already appeared in [9]. 

Theorem 1.2. The first nine coefficients of Ns{d) are given by 



3 6 



_ 5{5 - 1)(69(5^ - 85(5 + 92) ^2^.4 _ 5(8 - 1)((5 ^ 2)(702(5^ - 629(5 - 286) ^2^-5 
54 270 + 

(5((5 - 1)(^ - 2)(6028^^ - 15476(5^ + 11701(5 + 4425) ^ag-s 
+ 3240 + 

(5((5 - l)(g - 2) ((5 - 3) (13628(5^ - 6089(5^ - 295723 - 24485) 2^-7 
^ 11340 ^ 

(5(5 - 1)((5 - 2)((5 - 3)(282855(5'' - 931146(5=* + 4174905^ + 4252025 + 1141616)^25-8 

204120 ^ ' 

Let d*{6) denote the polynomiality threshold for Severi degrees, i.e., the smallest 
positive integer d* = d*{6) such that Ns{d) = N^'^ for d > d*. As mentioned above 
S. Fomin and G. Mikhalkin showed that d* < 26. We improve this as follows: 

Theorem 1.3. For 6 > I, we have d*{5) < S. 

In other words, A^"'''^ = Ns{d) provided d > 5 > I. L. Gottsche |8i Conjecture 4.1] 
conjectured that d* < [f] + 1 for S > 1. This was verified for 5 < 8 by S. Kleiman 
and R. Piene jH]. By direct computation we can push it further. 

Proposition 1.4. For 3 < 6 < 14, we have d*{6) = [f] + 1- 

That is, Gottsche's threshold is correct and sharp for 3 < 5 < 14. For 5 = 1,2 it 
is easy to see that d*{l) = 1 and d*{2) = 1. 

P. Di Francesco and C. Itzykson [6] hypothesized that d*{6) < 



(which is equivalent to 5 < — — — —). However, our computations show that this 
fails for 5 = 13 as d*{13) = 8. 

The main techniques of this paper are combinatorial. By the celebrated Correspon- 
dence Theorem of G. Mikhalkin fUi Theorem 1] one can replace the algebraic curve 
count by an enumeration of certain tropical curves. E. Brugalle and G. Mikhalkin 
[3111] introduced some purely combinatorial gadgets, called (marked) labeled floor di- 
agrams (see Section|2]), which, if counted correctly, are equinumerous to these tropical 
curves. Recently, S. Fomin and G. Mikhalkin [7] enhanced Brugalle and Mikhalkin's 
definition and introduced a template decomposition of labeled floor diagrams which 
is crucial in the proofs of all results in this paper, as is the reformulation of algebraic 
plane curve counts in terms of labeled floor diagrams (see Theorem 2.5). 



I + J26+, 
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This paper is organized as follows: In Section |2] we review labeled floor diagrams, 
their markings, and their relationship with the enumeration of plane algebraic curves. 
The proofs of Theorems and L2 are algorithmic in nature and involve a computer 



computation. We describe both algorithms in detail in Sections [3] and [5} respectively. 
The first algorithm computes the node polynomials Ns{d) for arbitrary 6, the second 
determines a prescribed number of leading terms of Ns{d). The latter algorithm 
relies on the polynomiality of solutions of certain polynomial difference equations: 
This polynomiality has been verified for pertinent values of 6 (see Section [sj. Propo- 



sition 



1.4 is proved in Section M by comparison of the numerical values of Ns{d) 



and iV^for various d and 6 (see Appendices |A] and |b]) . Theorem |l.3| is proved in 



Section H 

Competing Approaches: Floor Diagrams vs. Caporaso-Harris recursion. 

An alternative approach to computing the node polynomials Ns{d) combines poly- 
nomial interpolation with the Caporaso-Harris recursion [S]. Once a polynomiality 
threshold (io(5) has been established (i.e., once we have proved that Ns{d) = N'^'^ 
for d > do{6)), we can use the recursion to determine a sufficient number of Severi 
degrees N''-'^ for d > do{6), from which we then interpolate. 

This approach was first used by L. Gottsche [H Remark 4.1(1)]. He conjectured 
in Conjecture 4.1] the polynomiality threshold do{6) = [|] -|- 1, and combined it 
with the "Gottsche- Yau-Zaslow formula" [8| Conjecture 2.4] (now a theorem of Y.- 
J. Tzeng [12]) to calculate the putative node polynomials Ns{d) for 6 < 28. The 
Gottsche- Yau-Zaslow formula is a stronger version of polynomiality that allows one 
to compute each next node polynomial by calculating only two additional Severi 
degrees A^'^o{'5),5 g^j^^ j\jdoiS)+i,5 ^ which is done via the Caporaso-Harris formula. Since 
Gottsche's threshold dQ{S) = [|] -|- 1 remains open as of this writing, the algorithm 
he used to compute the node polynomials is still awaiting a rigorous justification. 

The first polynomiality threshold do{6) = 26 was established by S. Fomin and 
G. Mikhalkin [71 Theorem 5.1]. Using this result, one can compute Ns{d) for 6 <9 
but hardly any furthei]^ With the threshold (io(5) = S established in Theorem 1.3, 
it should be possible to compute Ns{d) for 5 < 16 or perhaps 6 < 17. 

By contrast, our Algorithm [T] does not involve interpolation nor does it require 
an a priori knowledge of a polynomiality threshold. Our computations verify the 
results of L. Gottsche's calculations for 6 < 14. In our implementations. Algorithm [T] 
is roughly as efficient as the interpolation method discussed above. (We repeat that 
the latter method depends on the threshold obtained using floor diagrams.) 

Gromov-Witten invariants. The Gromov-Witten invariant Nd^g enumerates irre- 
ducible plane curves of degree d and genus g through 3(i + (7 — 1 generic points in P^. 
Algorithm [1] (with minor adjustments, cf. Theorem 2.5[ 2)) can be used to directly 



compute Nd,g, without resorting to a recursion involving relative Gromov-Witten 
invariants a la Caporaso-Harris [5J 



We used an efE.cient C implementation of the Caporaso-Harris recursion by A. Gathmann. 
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Follow-up work. By extending ideas of S. Fomin and G. Mikhalkin [7J and of the 
present paper, we can obtain polynomiality results for relative Severi degrees, the 
degrees of generalized Severi varieties (see [5l[Tl]). This is discussed in the separate 



paper p]|; see Remark 3.9 



A. Gathmann, H. Markwig and the author |2] defined Psi-floor diagrams which 
enumerate plane curves satisfying point and tangency conditions as well as condi- 
tions given by Psi-classes. We prove a Caporaso-Harris type recursion for Psi-fioor 
diagrams, and show that relative descendant Gromov-Witten invariants equal their 
tropical counterparts. 

Acknowledgements. 1 am thankful to Sergey Fomin for suggesting this problem 
and fruitful guidance. 1 also thank the anonymous referee, Erwan Brugalle, Grigory 
Mikhalkin and Gregg Musiker for valuable comments and suggestions. Part of this 
work was accomplished at the MSRl (Mathematical Sciences Research Institute) in 
Berkeley, CA, USA, during the semester program on tropical geometry. I thank 
MSRI for hospitality. 

2. Labeled Floor Diagrams 

Labeled floor diagrams are combinatorial gadgets which, if counted correctly, enu- 
merate plane curves with certain prescribed properties. E. Brugalle and G. Mikhalkin 
introduced them in j3] (in slightly different notation) and studied them further in 
To keep this paper self-contained and to fix notation we review them and their 
markings following [7] where the framework that best suits our purposes was intro- 
duced. 

Definition 2.1. A labeled floor diagram V on a vertex set {1, . . . ,d} is a directed 
graph (possibly with multiple edges) with positive integer edge weights w{e) satisfy- 
ing: 

(1) The edge directions respect the order of the vertices, i.e., for each edge i — )■ j 
of V we have i < j. 

(2) (Divergence Condition) For each vertex j of V, we have 

def 



1 e edges e 

. e , . e . 

J -> A; « -> J 

This means that at every vertex of V the total weight of the outgoing edges is larger 
by at most 1 than the total weight of the incoming edges. 

The degree of a labeled floor diagram P is the number of its vertices. It is connected 
if its underlying graph is. Note that in [7] labeled floor diagrams are required to 
be connected. If V is connected its genus is the genus of the underlying graph 
(or the first Betti number of the underlying topological space). The cogenus of 
a connected labeled floor diagram V of degree d and genus g is given by 5{V) = 
{d-i)(d-2) _ ^ li T) is not connected, let di, ^2, • • • and Si, S2, ■ ■ ■ be the degrees 
and cogenera, respectively, of its connected components. Then the cogenus of P is 
'^jSj + "^j^ji djdj'. Via the correspondence between algebraic curves and labeled 
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floor diagrams (f3 Theorem 3.9]) these notions correspond hterally to the respective 
analogues for algebraic curves. Connectedness corresponds to irreducibility. Lastly, 
a labeled floor diagram V has multiplicit])^ 

m = n ^(^)'- 

edges e 

We draw labeled floor diagrams using the convention that vertices in increasing 
order are arranged left to right. Edge weights of 1 are omitted. 

Example 2.2. An example of a labeled floor diagram of degree d = 4, genus g = 1, 
cogenus 5 = 2, divergences 1, 1, 0, —2, and multiplicity /x = 4 is drawn below. 



O 



-O 



To enumerate algebraic curves via labeled floor diagrams we need the notion of 
markings of such diagrams. 

Definition 2.3. A marking of a labeled floor diagram V is deflned by the following 



three step process which we illustrate in the case of Example 2.2 



Step 1: For each vertex j of V create 1 — div{j) many new vertices and connect 
them to j with new edges directed away from j. 

O -O-^ 

Step 2: Subdivide each edge of the original labeled floor diagram V into two 
directed edges by introducing a new vertex for each edge. The new edges inherit 
their weights and orientations. Call the resulting graph T). 

2 _ 2 





Step 3: Linearly order the vertices of V extending the order of the vertices of the 
original labeled floor diagram V such that, as before, each edge is directed from a 
smaller vertex to a larger vertex. 



2 _ 2 




The extended graph V together with the linear order on its vertices is called a 
marked floor diagram, or a marking of the original labeled floor diagram P. 

We want to count marked floor diagrams up to equivalence. Two markings T>i, 
V2 of a labeled floor diagram V are equivalent if there exists an automorphism of 
weighted graphs which preserves the vertices of V and maps Vi to 1)2. The number 
of markings I'i'D) is the number of marked floor diagrams V up to equivalence. 



If floor diagrams are viewed as floor contractions of tropical plane curves this corresponds to 
the notion of multiplicity of tropical plane curves. 
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Example 2.4. The labeled floor diagram V of Example 2.2 has i^iT)) = 7 markings 
(up to equivalence): In step 3 the extra 1-valent vertex connected to the third white 
vertex from the left can be inserted in three ways between the third and fourth white 
vertex (up to equivalence) and in four ways right of the fourth white vertex (again 
up to equivalence). 

Now we can make precise how to compute Severi degrees N'^'^ and Gromov-Witten 
invariants N^^g in terms of combinatorics of labeled floor diagrams, thereby reformu- 
lating the initial question of this paper. Part |2] is not needed in the sequel and only 
included for completeness. It first appeared in [3l Theorem 1]. 

Theorem 2.5. [7, Corollary 1.9, Theorem 1.6] 

(1) The Severi degree N^'^ , i.e., the number of possibly reducible nodal curves in 
of degree d with 6 nodes through ^^^^^^ — § generic points, is equal to 



V 



where V runs over all possibly disconnected labeled floor diagrams of degree 
d and cogenus 5. 

(2) The Gromov-Witten invariant N^^g, i.e., the number of irreducible curves in 
of degree d and genus g through 3d + g — 1 generic points, is equal to 



V 

where V runs over all connected labeled floor diagrams of degree d and genus g. 



3. Computing Node Polynomials 

In this section we give an explicit algorithm that symbolically computes the node 
polynomials Ns{d), for given 6 > 1. (As A^'^''^ = 1 for > 1, we put NQ{d) = 1.) An 



m 



implementation of this algorithm was used to prove Theorem LI and Proposition 1.4 
We mostly follow the notation in [TJ Section 5]. First, we rephrase Theorem 1.1 
more compact notation. For 6 < 8 one recovers ^ Theorem 3.1]. For 5 < 14 this 
coincides with the conjectural formulas of [8| Remark 2.5]. 

Theorem 3.1. The node polynomials Ns{d), for 6 < 14, are given by the generating 
function "^syo ^^(d)x^ via the transformation 



J2Ns{d)x^ = exp lj2Qs{d) 

S>0 ^ S>1 



)X^ 
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where 

Qi{d) = 3{d-lf, 



Q2{d) = =f{d-l){Ud-25), 

Q3(d) = i (690^2 - 2364d + 1899), 

Q^(d) = i (-12060^2 + 47835d - 45207), 

Q5{d) = 1(217728(^2 - 965646(i + 1031823), 

Qg(d) = i(-4010328d2 + 19451628d - 22907925), 

Qj{d) = ^(74884932(i2 - 391230216d + 499072374), 

Qs(d) = |(-1412380980d2 + 7860785643d - 10727554959), 

Qsid) = i(26842726680ci2 - 157836614730d + 228307435911), 
Qioid) = (-513240952752^2 + 3167809665372d - 4822190211285), 
Qu{d) = 3^(9861407170992^2 _ 63560584231524(i + 101248067530602), 
Quid) = 3^ (-190244562607008^2 + 1275088266948600d - 2115732543025293), 
Qi3{d) = 3^(3682665360521280^2 - 25576895657724768d + 44039919476860362), 
Quid) = 3^ (-71494333556133600^2 + 513017995615177680d - 913759995239314452). 

In particular, all Qs{d), for 1 < 5 < 14, are quadratic in d. 

L. Gottsche conjectured that all Qsi^d) are quadratic. This theorem proves his 
conjecture for 5 < 14. 

The basic idea of the algorithm (see |7[ Section 5]) is to decompose labeled floor 
diagrams into smaller building blocks. These gadgets will be crucial in the proofs of 
all theorems in this paper. 

Definition 3.2. A template F is a directed graph (with possibly multiple edges) on 
vertices {0, . . . , /}, for / > 1, and edge weights w{e) G Z>o, satisfying: 

(1) If i — 7- j is an edge then i < j. 

(2) Every edge i A- i + 1 has weight w{e) > 2. (No "short edges.") 

(3) For each vertex j, 1 < j < / — 1, there is an edge "covering" it, i.e., there 
exists an edge z — )■ /c with i < j < k. 

Every template F comes with some numerical data associated with it. Its length 
l{r) is the number of vertices minus 1. The product of squares of the edge weights 
is its multiplicity /^(F). Its cogenus S{T) is 



(j - i)w{e) - 1 



For 1 < j < /(F) let Xj = Xj{T) denote the sum of the weights of edges i ^ k 
with i < j < k and define 

kmm(X) = max(xj - j + 1). 
i<j<i ■' 
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Sir) £(r) /i(r) ^(r) x(r) A;„,i„(r) p(r,A;) s{t) 




1 

2 



4 
1 




1 



(2) 



k-i 

2k + 1 



1 





Cr O O T) 



2 
2 
2 
2 
2 
2 
2 



1 
1 
2 
2 
2 
3 
3 



9 

16 

1 

4 

4 

1 

1 





1 
1 

1 
1 



(3) 
(4) 
(2,2) 

(3,1) 
(1,3) 



;i,i,i) 



'1,2,11 



3 
4 
2 
3 
2 
1 
1 



k-2 

(?) 

2k{k - 2) 
2k{k - 1) 
3(A; + 1) 
k{Ak + 5) 



1 
2 

1 
1 





Figure 1. The templates with 6(T) < 2. 



This makes k^i^i^) the smallest positive integer k such that F can appear in a floor 
diagram on {1, 2, ... } with left-most vertex k. Lastly, set 

1 if all edges arriving at I have weight 1, 
otherwise. 



^(r) 



Figure [T] (Figure 10 taken from [7]) lists all templates F with S(T) < 2. 

A labeled floor diagram V with d vertices decomposes into an ordered collection 
(Fi, . . . , Tm) of templates as follows: First, add an additional vertex d+1 {> d) to V 
along with, for every vertex j of P, l — div{j) new edges of weight 1 from j to the new 
vertex d + 1. The resulting floor diagram V has divergence 1 at every vertex coming 
from P. Now remove all short edges from V, that is, all edges of weight 1 between 
consecutive vertices. The result is an ordered collection of templates (Fi, . . . ,Fm), 
listed left to right, and it is not hard to see that ^^SiTi) = S{V). This process is 



reversible once we record the smallest vertex fcj of each template Fj (see Example 3.3 ) 



Example 3.3. An example of the decomposition of a labeled floor diagram into 
templates is illustrated below. Here, ki = 2 and ^2 = 4. 



V 



V 



O — M> 



o 



o — *-o- 




t 



(ri,r2)= (o) 



o— 
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To each template T we associate a polynomial that records the number of "mark- 
ings of F:" For k G Z>o let F(fc) denote the graph obtained from F by first adding 
k + i — l — Xi short edges connecting z — 1 to i, for 1 < i < /(F), and then subdividing 
each edge of the resulting graph by introducing one new vertex for each edge. By |7], 
Lemma 5.6] the number of linear extensions (up to equivalence) of the vertex poset 
of the graph F(fc) extending the vertex order of F is a polynomial in k, if k > k^[^{r), 
which we denote by P{T,k) (see Figure [T|. The number of markings of a labeled 
floor diagram V decomposing into templates (Fi, . . . , F^) is then 

m 

u{v) = l[Pir,,ki), 

1=1 

where ki is the smallest vertex of Fj in V. The algorithm is based on 

Theorem 3.4 ([7J, (5.13)). The Severi degree N'^^\ for d,6 > 1, is given by the 
template decomposition formula 

m d-l{r^)+e{rm) fc2-«(ri) 

(Pi , . . .^T'fji ) i 1 — ^min (Tm ) ^1 — ^min (Tl ) 

where the first sum is over all ordered collections of templates (Fi, . . . ,Tm), for all 
m > 1, with YlJLi^i^i) ~ ^> '^'^'^ sums indexed by ki, for 1 < i < m, are over 

^min(Fj) ^ ki ^ ki^i l(Ti), 



Expression (3.1) can be evaluated symbolically, using the following two lemmata. 
The first is Faulhaber's formula [10] from 1631 for discrete integration of polyno- 
mials. The second treats lower limits of iterated discrete integrals and its proof is 
straightforward. Here Bj denotes the jth Bernoulli number with the convention that 

Lemma 3.5 ([IQ]). Let f{k) = X]f=o'^*^* ^ polynomial in k. Then, for n > 0, 



(3.2) F(n) -Mj^nk) = E 7TT ^ ( ' D 

k=0 s=0 j=0 V / 

In particular, deg{F) = deg{f) + 1. 

Lemma 3.6. Let f{ki) and g{k2) be polynomials in ki and k2, respectively, and 
let ai,bi,a2,b2 E Z>o. Furthermore, let F{k2) = Ylkl=^ai f i^^) ^ discrete anti- 
derivative of f{ki), where k2 > ai + bi. Then, for n > max(ai + &i + b2, 02 + ^2); 

n—b2 k2—bi n— 62 

Yl 9{k2) Y /(^i) = E 9{k2)F{k2). 

k2=a2 fci=ai fc2=max(ai +61,02) 



Example 3.7. An illustration of Lemma 3.6 is the following iterated discrete inte- 
gral: 

fe2 = lfcl = l k2 = l ^ ■' k2=3 
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Data: The cogenus 6. 


Result: 


The node polynomial Ns{d). 


begin 






Generate all templates F with 6{r) < S; 




Ns{d) ^ 0; 




forall the ordered collections of templates T = (Fi, . . . , F^) with 




Er=i'^(r.) = 5do 






i 

Q 








while i < m do 








ai ^ 








max (A;min(Fi), A;min(Fi_i) +/(Fj_i), . . . , /Cmin(Fi) + /(Fi) H h/(Fi_i)); 






end 






while i < m — 1 do 
















i i- i + 1; 






end 


















Nsid) ^ Nsid) + Q^d); 




end 




end 







Algorithm 1: Algorithm to compute node polynomials. 



Using these results Algorithm [T] computes node polynomials Ns{d) for an arbitrary 
number of nodes S. The first step, the template generation, is explained later in this 
section. 

Proof of Correctness of Algorithm^ The algorithm is a direct implementation of 
Theorem |3.4[ The m-fold discrete integral is evaluated symbolically, one sum at a 



time, using Faulhaber's formula (Lemma 3.5). The lower limit of the ith. sum is 
given by an iterated application of Lemma 3.6[ □ 



As Algorithm [T] is stated its termination in reasonable time is hopeless for 6 > 8 
or 9. The novelty of this section, together with an explicit formulation, is how to 
implement the algorithm efficiently. This is explained in Remark |3.8[ 



Remark 3.8. The running time of the algorithm can be improved vastly as follows: 
As the limits of summation in (3.1) only depend on k^miXi), K^i) ^i^m), we 
can replace the template polynomials PiVi^ki) by ^P(Fj,fcj), where the sum is 
over all templates Fj with prescribed (fcmin,^^^)- After this transformation the first 
sum in (3.1) is over all combinations of those tuples. This reduces the computation 



drastically as, for example, the 167885753 templates of cogenus 14 make up only 



343 equivalence classes. Also, in (3.1) we can distribute the template multiplicities 
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fj.{Ti) and replace P{Ti,ki) by fi{Ti)P(Ti, ki) and thereby eliminate n/^(ri)- An- 
other speed-up is to compute all discrete integrals of monomials using Lemma |3.5 
advance. 



m 



The generation of the templates is the bottleneck of the algorithm. Their number 
grows rapidly with 6 as can be seen from Figure [3j However, their generation can 
be parallelized easily (see below). 

Algorithm [l] has been implemented in Maple. Computing Nii(d) on a machine 
with two quad-core Intel(R) Xeon(R) CPU L5420 @ 2.50GHz, 6144 KB cache, and 
24 GB RAM took about 70 days. 

Remark 3.9. Using the combinatorial framework of floor diagrams one can show 
that also relative Severi degrees (i.e., the degrees of generalized Severi varieties, see [S, 
E]) are polynomial and given by "relative node polynomials" [D, Theorem 1.1]. This 
suggests the existence of a generalization of Gottsche's Conjecture [H Conjecture 2.1] 
and the Gottsche-Yau-Zaslow formula [8j Conjecture 2.1]. Thus, the combinatorics 
of floor diagrams lead to new conjectures although the techniques and results seem 
to be out of reach at this time. 

Remark 3.10. We can use Algorithm [l] to compute the values of the Severi degrees 
jsjd,5 prescribed values of d and 5. After we specify a degree d and a number of 
nodes 5 all sums in our algorithm become finite and can be evaluated numerically. 
See Appendix [b| for all values of A^'^-'' for < 5 < 14 and 1 < < 13. 



Proof of Proposition l.^^. For 1 < 5 < 14 we observe, using the data in Appendices 
A|and|B| that N5{d) = N'^'^ for all do{5) < d < S, where do{5) = [f] +1 is Gottsche's 
threshold. Furthermore, Ns{do{S) - 1) ^ N'^o(S)-i,s all 3 < 5 < 14. □ 



Template Generation. To compute a list of all templates of a given cogenus one 
can proceed as follows. First, we need some terminology and notation. An edge 
i — J- j of a template is said to have length j — i. A template F is of type a = (ajj), 
i,j G Z>o, if r has aij edges of length i and weight j. Every type a satisfies, by 
definition of cogenus of a template, 

(3.3) 5^«.,(^-j-l)=<5(r). 

Note that an = as short edges are not allowed in templates. The number of types 
constituting a given cogenus 6 is finite. 

We can generate all templates of type a using a branch-and-bound algorithm which 
slides edges in a suitable order. Let Fq be the unique template of type a with all 
edges emerging from vertex 0. Fix a linear order on the set of edges of type a. For 
[0 ll 



example, if a 



2 



we could choose: 



o — *~o 

2 
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Data: A graph A with a distinguished edge Ci. 
Result: An infinite directed tree of graphs with root A. 
begin 

forall the edges 62 of A with 62 > ei (in the fixed order) do 
B ^ graph obtained from A by moving 62 to the next vertex; 
if the natural partial order (from left to right) of the edges of B that are 
of the same type as 62 is compatible with the fixed order then 
Insert B child of A; 
Execute this procedure with input (5,62); 
end 
end 
end 

Algorithm 2: A recursion which can generates a tree containing all templates 
of a given type. 



Algorithm [2] applied to the pair (Fq, cq), where cq is the smallest edge of Fq, creates 
an infinite directed tree with root Fq all of whose vertices correspond to different 
graphs. Eliminate a branch if either 

(1) no edge of the root of the branch starts at vertex 1, 



or 



(2) condition (3) in Definition 3.2 is impossible to satisfy for graphs further down 
the tree. 

1" 
2 



See Figure 



for an illustration for a 




Figure 2. Branch-and-bound tree for a 



1 
2 



COMPUTING NODE POLYNOMIALS FOR PLANE CURVES 



13 



A complete, non-redundant list of all templates of type a is then given by all 



remaining graphs which satisfy condition (3) of Definition 3.2 as every template can 
be obtained in a unique way from Fq by shifting edges in an order that is compatible 
with the order fixed earlier. Note that it can happen that a non-template graph 
precedes a template within a branch. For an example see the graph in brackets in 
Figure |3} Template generation for different types can be executed in parallel. The 
number of templates, for 5 < 14, is given in Figure [3] 



b 


of templates 


b 


of templates 


b 


of templates 


1 


2 


6 


1711 


11 


2233572 


2 


7 


7 


7135 


12 


9423100 


3 


26 


8 


29913 


13 


39769731 


4 


102 


9 


125775 


14 


167885753 


5 


414 


10 


529755 







Figure 3. The number of templates with cogenus b < 14. 



4. Threshold Values 

S. Fomin and G. Mikhalkin [71 Theorem 5.1] proved polynomiality of the Severi 
degrees iV'^''^ in d, for fixed b, provided d is sufficiently large. More precisely, they 
showed that Ns{d) = N"^'^ for d > 2b. In this section we show that their threshold 



can be improved to d > b (Theorem 1.3). 



We need the following elementary observation about robustness of discrete anti- 
derivatives of polynomials whose continuous counterpart is the well known fact that 
Ja~i ^ f{x)dx = if f{x) = on the interval (a — s — 1, a — 1). 

Lemma 4.1. For a polynomial f{k) and a G Z>o let F{n) = Yl^=a /(^) poly- 
nomial in n uniquely determined by large enough values of n. (F{n) is a polynomial 



by Lemma 3.5 ) If we have /(a— 1) = ■ ■ ■ = f(a—s) = for some < s < a (this con- 
dition is vacuous for s = 0) then it also holds that F{a — 1) = ■ ■ ■ = F{a — s — 1) = 0. 
In particular, X]fc=a /(^) ^■^ ^ polynomial in n, for n > a — s — 1. 

Even for s = the lemma is non-trivial as, in general, F{a — 2) 7^ 0. 

Proof. Let G{n) be the polynomial in n defined via G{n) = X]fc=o /(^) large n. 
Then F{n) = G{n) — ^^Iq f{k) for all n G Z>o. In particular, for any < z < s, we 
have F{a-i-l) = G{a-i-l)- YX=l /(^) = G{a~i-l)- /(^) = 0- □ 

Recall that for a template F, we defined fcmm = fcmin(r) to be the smallest A; > 1 
such that k + j — 1 > Xj(F) for all 1 < j < ^(F). Let jo be the smallest j for which 
equality is attained (it is easy to see that equality is attained for some j). Define 
s(F) to be the number of edges of F from jo — 1 to jo (of any weight). See Figure [ijfor 
some examples. The following lemma shows that the template polynomials P(F, k) 



satisfy the condition of Lemma 4.1 
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Lemma 4.2. With the notation from above it holds that 

p(r, k^,^ - 1) = p(r, kmin - 2) = ■ ■ ■ = p(r, k^,^ - s{t)) = o. 



Proof. Recall from Section|3]that, for k > /cjnin(r), the polynomial P{T, k) records the 
number of linear extension (up to equivalence) of some poset r(fc) which is obtained 
from r by first adding k + j — 1 — >Cj{T) "short edges" connecting j — 1 to j, for 
1 < j < /(r), and then subdividing each edge of the resulting graph by introducing 
a new vertex for each edge. 

Using the notation from the last paragraph notice that /cmin + jo — 1 = ^io(r), and 
thus r(fc) has k — fcmin "short edges" between jo — 1 and jo- Every linear extension of 
r(fc) can be obtained by first linearly ordering the midpoints of these k — k^i,^ "short 
edges" and the midpoints of the s{T) many edges of V connecting jo — 1 and jo 
before completing the linear order to all vertices of '^{k)- Therefore, the polynomial 
{k — /cjnin + 1) ■ ■ ■ (/c — /Cmin + ^(r)) dividcs P{T, k) . □ 



Before we can prove Theorem 1.3 we need a last technical lemma. 

Lemma 4.3. Using the notation from above we have, for each template T , 

k^^{T) - s{T) + 1{T) - e{T) < 5(1) + 1. 

Proof. As before, let jo be the smallest j in {1, . . . , l(T)} with kmin + j — 1 = >'^j(X)- 
It suffices to show that Xjo(T) — jo — s(r) + l(T) — e(T) < S(T). 

Let r' be the template obtained from F by removing all edges i k with either 
k < joOY i> jo. It is easy to see that 1{T) - e{T) - (/(!') - e{T')) < 6(T) - 6{T'). 
Thus, we can assume without loss of generality that all edges i — )■ of F satisfy 
i < io < k. Therefore, as XjoiV) = ^ ^^^^^ ^ of r ^^(^) suffices to show that 

(4.1) /(r)-£(r)< Yl [wt(e)(len(e)-l)-l] +s(r)+jo, 

edges e of r 

where len(e) is the length of an edge e. The contribution of the s(r) edges of F 



between jo — 1 and jo to the sum is — s(r), thus the right-hand-side of (4.1) equals 
(4.2) ^[wt(e)(len(e)-l)-l] +jo 

with the sum now running over all edges of F of length at least 2. If there are no 
such edges, then /(F) = 1 and we are done. Otherwise, if e(T) = 1, expression ( 4.2[ ) 



equals ^{\en{e) - 2) + jo, which is > /(F) - 2 + jo or > /(F) - 3 + jo if jo e {1, /(F)} 
or 1 < jo < /(F), respectively (by considering only edges adjacent to vertices and 
/(F) of F). In either case the result follows. 



If e{T) = then expression (4.2) is > /(F) + (/(F) - 3 + jo) or > /(F) - 2 + j o if 



jo € {1,/(F)} or 1 < jo < /(F), respectively. This completes the proof. □ 
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Proof of Theorem 1.3. By Lemma 3.6 and repeated application of Lemmata 4.1 



and |4.2| it suffices to show that d> 6 simultaneously implies 

d ^ /(Lf^) ^(r^) -|- fcmin(rm) ^{,^m) 1; 

~l~ ^min (r?Ti— 1 ) ^ij^m—l] 

(4.3) 



d > /(r„) - £(r„) + /(r„_i) + • • • + /(Li) + A;^in(ri) - s(ri 

for all collections of templates (Fi, . . . , Tm) with YllLi ^O^i) = ^■ 

The ffist inequality is a direct consequence of Lemma |4.3[ For the other inequali- 
ties, notice that Z(F) — e{r) < 5(F) for all templates F, hence 

/(FJ-£(FJ-1<5(FJ-1 

and 

l(T,) -1 <S(Ti), foT2<i<m-l. 



By Lemma 4^ we have 

/(Fi) + k^UTi) - s(Fi) - 1 < 5(Fi) + 1 



as e{Ti) < 1, and the right-hand-side of the last inequality of (4.3) is < Xli^i ^(^i) = 
6 < d. The proof of the other inequalities is very similar. □ 

5. Coefficients of Node Polynomials 
The goal of this section is to present an algorithm for the computation of the 



coefficients of Ns(d), for general 6. The algorithm can be used to prove Theorem L2 



and thereby con&m and extend a conjecture of P. Di Francesco and C. Itzykson in 
|6] where they conjectured the 7 terms of Ns{d) of largest degree. 

Our algorithm should be able to find formulas for arbitrarily many coefficients of 
Ns{d). We prove correctness of our algorithm in this section. The algorithm rests on 



the polynomiality of solutions of certain polynomial difference equations (see (5.7)). 



First, we fix some notation building on terminology of Section [3] By Remark |3.8 
we can replace the polynomials P(F, k) in (3.1) by the product /i(F)P(F, k), thereby 
removing the product Yl f^i^i) of the template multiplicities. In this section we write 
P*(F, k) for /i(F)P(F, k). For integers i > and a > let Mj(a) denote the matrix 
of the linear map 

n-l{r) 

(5.1) f{k)^ Yl E P*i^,k)-f{k), 

r:<5(r)=ifc=A;n,in(r) 

where f{k) = Cok"" + Cik"""^ + ■ ■ ■ , a polynomial of degree a, is mapped to the 



polynomial Mi{a){f{k)) = don'^'^^'^^ + rfin^^* + ■ ■ ■ in n. (By Lemma 3.5 and the 



proof of Lemma 5.1 the image has degree a + i + 1. 
define Mf^ia 



Hence Mj(a)c = d. Similarly, 



(5.2) 



to be the matrix of the linear map 

n-/(r)+e{r) 

/(^)^ E E P*i^,k).f{k). 

r:<s(r)=i fc=fc^in(r) 
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Later we will consider square sub-matrices of Mj(a) and M, 



end/ 



by restriction to 



the first few rows and columns which will be denoted Mj(a) and Mf (a) as well. 
Note that Mj(a) and Mf'^'^{a) are lower triangular. For example, for a large enough, 



Ml (a) 



40 



6 

a+2 
5a+8 
a+1 

|a + 3 
l{Aa + l)a 





6 
a+1 
5a+3 

2O + 2 



;i3a2 - 20a + 7)a -a^ + |a - | 






6 

5a- 2 
' a-1 









a— 1 

5a-7 6 



Lemma 5.1. The first a + i rows of Mi{a) and M/"'^(a) are independent of the lower 



limits of summation in (5.1) and (5.2), respectively 



Proof. It is an easy consequence of the proof of [71 Lemma 5.7] that the polynomial 
P*{T,k) associated with a template T has degree < S(T). Equality is attained by 
the template T on vertices 0, 1, 2 with i edges connecting and 2 (so 6(T) = i). As 
discrete integration of a polynomial increases the degree by 1 the polynomial on the 
right-hand-side of (5.1) has degree 1 + i + a. □ 



The basic idea of the algorithm is that templates with higher cogenera do not 
contribute to higher degree terms of the node polynomial. With this in mind we 
define, for each finite collection (Fi, . . . ,Fm) of templates, its type r = (r2,r3, . . . ), 
where Tj is the number of templates in (Fi, . . . , F^) with cogenus equal to i, for i >2. 
Note that we do not record the number of templates with cogenus equal to 1. 

To collect the contributions of all collections of templates with a given type r, 
let T = (r2,r3,...) and fix 5 > Ylij>2'^3 (^^ ^^^^ there exist template collections 
(Fi,...,Fm) of type r with — We define two (column) vectors Cr{5) 

and C'^^'^{5) as the coefficient vectors, listed in decreasing order, of the polynomials 

(5.3) Yl E P*i'^m,km)--- ^*(ri,fcl) 

km — ^min (r™) ki=k 

min (ri) 

and 

n-«(r„) + £(r) krr^-liTm-l) ki-liVi) 

(Ti,. ■■1^771.) km—ki-niYi{T^rn) ^m — 1 — ^miii (^771 — 1 ) ^1— ^min(ri) 

in the indeterminate n, where the respective first sums are over all ordered collections 
of templates of type r. 
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It might look like Cr{S) is a product of some matrices Mj(a) applied to the poly- 
nomial 1. However, this is not the case. For example, note that 



(0,0,...) 



(2) 



r 9 -| 




r 9 -1 


2 




2 


-34 




-34 


88 

179 




88 

179 


2 


2 


30 




27 











Ml (2) ■ Mi(0) 



This is because, when iterated discrete integrals are evaluated symbolically, the lower 
limits of integration of the outer sums can change depending on the limits of the inner 



sums (cf. Lemma 3.6). This observation makes it necessary to compute initial values 
for recursions (described later) up to a large enough 6. 

Before we can state the main recursion we need two more notations. For a type 
r = (r2, r3, . . . ) and i > 2 with Tj > define a new type rj,j via (rj,j)j = Tj — 1 and 
(rj,j)j = Tj for j 7^ i. Furthermore, let def(r) = J2j>2U ~ be the defect of r. 
The following lemma justifies this terminology. 



Lemma 5.2. The polynomials (5.5) and (5.4) are of degree 26 — def(r). 



Proof. Let (Fi, . . . , F^) be a collection of templates with X^i ^(Xi) = ^ tyP^ t. 
Then, by applying the argument in the proof of Lemma [5. 1| to each Fj, the polyno- 



mials (5.3) and (5.4) have degree S + m. The result follows as 



def(r) = 5^5(F,)- 




i=l 






m 


E 


[(^ 


i=l 


i>2 


_ \i:<5(r,) = 







m. 



J>2 



□ 



The last lemma makes precise which collections of templates contribute to which 
coefficients of Ns{d). Namely, the first coefficients of Ns{d) of largest degree 
depend only on collections of templates with types r such that def(r) < A^. The 
following recursion is the heart of the algorithm. 

Proposition 5.3. For every type t and integer 6 large enough, it holds that 
Cr{6) = ^ M,(25 - z - 1 - defir))Crd6 - t) 

(5.5) i-.Tiy^Q 

+ Mi{26 - 2 - dei{T))Cr{6 - 1). 
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More precisely, if we restrict all matrices Mi to be square of size N — def (r) and all 



Cr to he vectors of length N — def (r), then recursion (5.5) holds for 



5 > max 



N + i 




Proof. The coefficient vector Cr{S) is defined by a sum that runs over all collections 
of templates (Fi, . . . , r^) of type r (see (5.3)). Partition the set of such collections 
by putting 6(Tm) = 1, or 6(Tm) = 2, and so forth. This partitioning splits expression 



(5.3) exactly as in (5.5). 



A summand can be written as a product of some matrix Mj and some vector C-j-ii 
if 6 is large enough, namely if Mj does not depend on the lower limits in (5.3). If we 



can factor then the polynomials (5.3) defining Crii{S — i) and Cr{5 — 1) have degrees 
2(5 -i)- def (4i) = 25 - 2i - def(r) + (i - 1) = 25 - i - 1 - def(r) 



by Lemma 5.2 and, similarly, 25 — 2 — def(r), respectively. By Lemma 5.1, if the 
matrix Mj(25 — i — 1 — def(r)) is of size N — def(r), then it does not depend on the 
lower limits if and only if 5 > In order for Cr[S) to be defined (and the above 
identity to be meaningful) we need to impose 5 > Ylj>2 3'^j- '— ' 

Remark 5.4. Later, when we formulate the algorithm, we need to solve recursion 



(5.5) together with an initial condition in order to obtain an explicit formula for the 



first N — def (r) entries of C^(5). It suffices to take 



(5.6) 



5o(r) =^ max 



'N 




as for any 5 > 5o(r) the vector Cr{S) of length N — def(r) can be written in terms 
of matrices Mi and vectors Cr>{S') for various types r' and integers 5' < 5. 

We propose Algorithm |3] for the computation of the coefficients of the node poly- 



nomial Ns{d). We explain how to solve recursion (5.5) below 



Proof of Correctness of Algorithm^ Proposition 5.3 guarantees that (7^(5) is uniquely 
dete rmin ed by recursion (5.3). By a similar argument as in the proof of Prop osi- 
tion 5.3 we see that C^^'^{5) is given by the formula in Algorithm p| By Lemma 5.2 



all contributions of template collections of type r to the node polynomial Ns{d) are 
in degree 25 — def(r) or less. Hence, after shifting C^°'^(5) by def(r), their sum is 
the coefficient vector of Ns{d). □ 



To solve recursion (5.5) for a type r we make use of the following (conjectural) 
structure about Cr{6) which has been verified for all types r with def(r) < 8. This 
refines an observation of L. Gottsche |H1 Remark 4.2 (2)] about the first 28 (conjec- 
tural) coefficients of the node polynomial Ns{d). 



Conjecture 5.5. All entries of Ct-{5) are of the form ^ times a polynomial in 5. 
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Data: A positive integer N. 

Result: The coefficient vector C of the ffist coefficients of Ns{d). 
begin 

Compute all templates T with 6(T) < N; 
forall the types t with def(r) < do 



Compute initial values Ct-((5o(t)) using (5.3), with So{t) as in (5.6); 



Solve recursion ( |5.5[ ) for ffist A^ — def(r) coordinates of Cr{S); 
Set 

C-d(5) ^ J2 Mr''{25 - z - 1 - def(r))a4.(5 - t) 

+ Mi^"<^(2(5 - 2 - def(r))a(5 - 1); 
end 

forall the types r with def(r) < A^ do 

Shift the entries of C^^'^{5) down by def(r); 
C shifted C^'^'i (5); 

end 



end 



Algorithm 3: Computation of the leading coefficients of the node polynomial. 



Now, to solve recursion (5.5), we ffist extend the natural partial order on the types 
r given by |r| = ^j>2 T?' a linear order with smallest element r = (0, 0, . . . ). For 
example, for A^ = 4, we could take 

(0,0,0) < (1,0,0) < (0,1,0) < (0,0,1) < (1,1,0) < (2,0,0) < (3,0,0). 

Then solve recursion ( 5.5[ ) for each r, in increasing order, using the lowertriangularity 



of the matrices Mj. For example, to compute the second entry 4fp((5) of Ci^i{6) 



(assuming Conjecture 5.5[ ), where p{5) is a polynomial in 5, we need to solve 
Ci,i((5) = Ml (25 - 5)Ci,i(5 - 1) + M2(25 - 6)Co,i(5 - 2) + M^{26 - 7)Ci,o(5 - 3), 
or, explicitly. 







"* " 








* * 





(<5-l) 



■p{5 - 1) 



+ 



"* " 








"* " 






* * 




* 


+ 


* * 




* 



The *-entries in the vectors Co,i and Ci^ are known by a previous computation. The 
*-entries in Mi, M2 and M3 are given by (5.3). The proof of Lemma 5.1 implies that 
all denominators of Mi{a) in row j are a+i—j+2 or 1 (after cancellation). To compute 
p(5), or, more generally, the jth entry in Ct-(5), we ffist clear all denominators and 
then solve the polynomial difference equation with initial conditions 

(25 - def (r) - j + l)3p(5) = p{5 - 1) + g(5). 
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where q{6) is a rather comphcated polynomial depending on earlier calculations and 



So{t) is as in (5.6). One way to solve (5.7) is to bound the degree of the polynomial 



p{6) and solve the corresponding linear system. 



Note that a difference equation of the form (5.7) need not have a polynomial 



solution in general. Conjecture 5.5 is equivalent to all recursions (5.7) appearing in 
Algorithm |3] to have a polynomial solution. 

As in Section [3| (Remark |3.8[ ), Algorithm [3] can be improved significantly by sum- 
ming the template polynomials P(T, k) for templates T with fixed [k^miX), ^(r), ^(r)) 
in advance. Algorithm [3] has been implemented in Maple. Once the templates are 
known the bottleneck of the algorithm is the initial value computation which, with an 
improved implementation, should be faster than the template enumeration. Hence 



we expect Algorithm |3] to compute the first 14 terms of Ns{d) in reasonable time. 

Appendix A. Node Polynomials for 5 < 14 
An explicit list of Ns{d), for 6 < 14, is as below. These polynomials are given 



implicitly in Theorem 3.1 For 6 < 8 this agrees with [9l Theorem 3.1]. For 6 < 14 



this coincides with the conjectural (implicit) formulas of [H Remark 2.5]. 

No{d) = 1, 

Ni{d) = 3{d - if , 

N2{d) = ^(d- 2){3d^ - 3d - 11), 

N3{d) = -d« - 27<i5 + ^d^ + l?5rf3 _ 229^2 _ ^^ + 525 
^ ^ 2 2 2 2 

N4d) = - 27d^ +^d^- 642^4 _ 252M3 + ?^d^ + ^^d - 8865, 

^ ^ 8 4 8 4 

81 ,1(1 81 ,„ 27 ,o 2349 . 127071 128859 59097 . 3528381 ,o 

Nsid) = —d^" d^ d^ H d^ - 1044d^ d^ + H d^ d^ 

40 4 8 4 20 8 2 40 

946929 

d+ 153513, 

20 

^,(,) = 51,- _ 243^,, _ 81^,0 + 8667^« _ 9297^3 _ 47727^, 2458629^, 3243249^, 

80 20 20 16 8 5 80 40 

6577679 . 25387481 6352577 ,0 8290623 , 

d^ H d^ H d - 2699706, 

20 80 4 20 

^,(,) = 243^,, _ 243^,3 _ 243^,, 30861^,, _ 38853^,o _ 802143^, 3140127^3 + 18650493^, 

560 40 80 80 40 80 80 140 

54903831 ,fi 72723369 124680069 . 213537633 ,0 3949576431 ,2 188754021 , 

d" d^ H d* H d-^ d^ d 

80 80 20 80 140 140 

+ 48016791, 

,w A 729 729,,. 243,-14 35721,-,., 25839,-,, 320841,,, 11847087,,,, 

Ns(d) = d^^ d^^ d^^ H d" d^^ d^^ H d^° 

4480 280 140 160 40 40 320 

170823033 ,„ 6685218 ,0 1758652263 ,7 1102682031 .r 59797545 ,^ 510928080111 

H d-^ d^ H d^ H d^ d^ 

1120 7 1120 80 8 4480 

3283674393 558215113803 3722027733 , 

d^ H d^ d- 861732459, 

1120 1120 56 

243 ,0 2187 ,7 729 ,fi 121743 99549 824823 ,0 8776593 ,„ 74122857 ,, 

Ng(d) = d^** d" d^" H d^^ d^* d" H d^'' H d" 

^ ' 4480 2240 896 1120 280 160 320 560 

2188424421 132610923 ,q 11404136871 ,« 2852923401 ,7 3523392270287 ,« 

d^" d^ H d** H d^ d" 

2240 70 560 224 13440 

4109675615 261844582229 ,4 2156232149611 29528525065861 ,, 438722045999 , 

H d^ H d^ d"' d^ H d 

448 128 3360 3360 168 

+ 15580950065, 
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' 44800 2240 2240 8960 4480 700 640 

205722099 4375229931 ,,o 38815692777 30958937073 3413568339 ,„ 

2240 5600 22400 1400 224 

_ 3624162885799 134470136581^7 27023302169081^6 22514488581251 811909836973903^4 
8960 2800 ^ 5600 8960 22400 

253124357071961 , 867510616107447 , 2800250331071 

H (P + (f d - 283516631436, 

11200 5600 40 

2187 2187 o, 729 ,„ 150903 303993 56670273 ,7 47717667 

Nii(d) = d-'-' d-^^ d^" H d^^ d^^ d" H d^" 

^ ' 492800 22400 6400 8960 4480 44800 5600 

295979589 ,5 11410430877 ,4 4051907631 ,3 52491198663 .j, 3418059518271 ,„ 

H d d d H d H d 

5600 22400 3200 2800 246400 

20587006282467 2236636275459 ,q 49175916627959 ,0 1464110674563 - 

44800 22400 6400 256 

1946239824069277 . 3767687640687823 ,5 14264414890838423 ,4 940418544772283 ,3 

d H o H d d 

22400 44800 22400 1600 

_ 168280746183263029^, + 5187507215325, 

61600 3080 

= ^H^,- _ ^,.3 _ 218^,.. ^ 496449^,, _ 136809^,„ _ 1618623^,, 674946837^,3 

^ ' 1971200 246400 61600 89600 5600 3200 179200 

2321658693 893195181 .g 34334301951 ,5 289702847403 ,4 1245724147341 ,3 

H a d d H o H d 

89600 3200 44800 22400 123200 

803786361621603 65497548165237 16192295343681 ,in 792669234543351 ,„ 

d H d H d d 

1971200 492800 1792 89600 

9506773589164709 ,0 6296062244021929 ,7 11029935159768347 ,r 582428855393100577 

H d'^ H dr d^ 

67200 33600 7168 268800 

5477484616918678589 10067756533588172119 ,., 4454424013895459501 ,2 

d* H d-^ H d^ 

492800 739200 92400 

111952943233924509 

d - 95376705265437, 

3080 

6561 6561 19683 1620567 ,0 88209 „ 3212703 262066023 

iVi3(d) = d^" d^-' d^^ H d^-^ d^^ d^^ H d^" 

^ ' 25625600 985600 1971200 985600 11200 17920 179200 

494726373 ,„ 673360047 35350103511 ,7 20952637821 ,f, 3013479294723 

-I cjJ^'^ d d H d H d 

44800 5120 89600 2800 492800 

_ 580214902388013^14 1666286215401123^13 16384163286402207^,2 _ 909876952033137 

1971200 12812800 1971200 89600 

7649416285706767 25855007471662161 ,q 65085797443981191 ,0 108443195356282427 ,7 

_ d^^ -I d H d — d 

44800 89600 25600 22400 

52991400162927629917 1976324604711031517 ,^ 13580753080243105219 . 

dy H d^ + d* 

1971200 39424 70400 

73274705967431063281 , 68173290776099374391 2813974748454890667 , 

d^ d^ H d+ 1761130218801033, 

246400 80080 3640 

N (d) = ^^^^3 o;28 _ _L9683_^27 _ 6561 1751787 ^2,, _ 4529277 ^24 _ 562059 ^23 

358758400 12812800 2562560 3942400 1971200 9856 

398785599^22 5214288411^2, _ 4860008991^20 _ 63174295089^,9 332872084467^,8 

788480 1254400 89600 358400 89600 

3103879378581^,7 _ 4913807521304691 899178800016807^,5 279086438050359453^,4 

^ 985600 27596800 8968960 44844800 

468967272863997483^,3 318443311640108577^,2 328351365725506869^,, 

51251200 1971200 985600 

112058081 i()8()r)71l)23 m 9448S()1()2S 1188 13949 „ 3088078r)21()7:-i() i()(ill3 « 

J d d d 

358400 1254400 689920 

444525313669622586903 .7 11429038221675466251 269709254062572016617 

H d^ H d'^ d'^ 

3942400 24640 246400 

74660630664748878665353 . 140531359469510983018159 , 16863931195154225977601 , 

d* H d-^ H d^ 

22422400 22422400 1121120 

64314454486825349085 

d - 32644422296329680. 

4004 
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Appendix B. Small Severi degrees 



Below we list the Severi degrees 

j^d,5 for < (5 < 14 and 1 < < 13, which were 
obtained by Algorithm [l] (also see Remark 3. 10 ) . Together with the node polynomials 
of Appendix [A] this is a full description of all Severi degrees N"^'^ for 6 < 14 



see 



Theorem 1.3 



The solid line segments indicate the polynomial threshol d d* {6) of 

The 



N'^'^. The dashed line segments illustrate the threshold of our Theorem 1.3 
Severi degrees A^'^''' in italic agree with the Gromov-Witten invariants N , (d-i){d-2) 

as ioT d > 5 + 2, every plane degree d curve with 5 nodes is irreducible. 
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